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> 1 ABSTRACT 

. Aims. We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. 
\^ • Methods. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 
£Nj ' 14.2 hours of good quality very high energy (VHE) y-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the 

* , Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy y-ray (HE, 0. 1 GeV-100 GeV) data from the Fermi 

satellite. 

Results. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE y-rays at a statistical significance of 5.6cr. A soft 
VHE spectrum with a photon index of -3.2 ± 0.6 was determined. No significant short term flux variations were found. We model the spectral 
energy distribution using a one-zone SSC-model, which can successfully describe our data. 

, , Key words. BL Lac objects: individual(B3 2247+381)-galaxies:active-gamma rays 
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1 . Introduction al. 2007a), BL Lacertae (part of the data taken as ToO due to op- 
tical high state, Albert et al. 2007b), 3C 279 (Albert et al. 2008a, 
The number of known extragalactic very high energy (VHE, E Aleksic et al. 201 la), MAGIC J0223+430 (Aliu et al. 2009), S5 
> 100 GeV) y-ray emitters has increased from 6 to almost 50 0716+714 (Anderhub et al. 2009), B3 2247+381 (Mariotti et 
(July 201 1) during the past five yearfU Most of these objects are al 2010i this paper) and most rece ntly 1ES 1215+303 (Mariotti 
active galactic nuclei (AGN), especially blazars (flat spectrum et al 2 011). This suggests that high optical states are an in- 
radio quasars and BL Lac objects), which are known for their dication of a h i gher sta te also in VHE y-ray band, at least in 
large variability across the electromagnetic spectrum from radio some sources . Additionally for PG 1553+1 13 a long term study 
to VHE y-rays. Eight of these new discoveries were made dur- (Aleksic, J., Alvarez, E.A., Antonelli, L.A., Antoranz, P. et al., 
ing high optical states of these sources, reported by the Tuorla A pJ, in prep.) suggests a correlation between the optical and 
blazar monitoring program, which have triggered MAGIC obser- VHE y-ray flux, while for PKS 2155-304 it seems that, at least, 
vations: Mrk 180 (Albert et al. 2006), 1ES 101 1+496 (Albert et in some cases the two wave bands are correlated (Foschini et al. 
2007; Aharonian et al. 2009). 



http://www.mpp.de/ rwagner/sources/ 
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The object B3 2247+381 was selected from the sample pre- 
sented in Nieppola et al. (2006) with a reported X-ray flux F(> 
IkeV) > 2mJy. Donato et al. (2001) classify it as high-energy- 
peaked BL Lac object and the reported X-ray flux of the source 
is F(>lkeV)= 0.6mJ>B B3 2247+381 (z = 0.119, Falco et al. 
1988) has been previously observed by the MAGIC telescope as 
part of the systematic search of VHE T-rays from X-ray bright 
BL Lac objects (Aleksic et al. 2011bfl MAGIC observed the 
source between August and September 2006, which resulted in 
an upper limit F(>140GeV) < 1.6 ■ 10 _1I cm _2 s _1 . The source 
has been monitored in the R-band by the Tuorla blazar monitor- 
ing program ever since. 

B3 2247+381 was also included in the list of potentially 
interesting TeV sources released to the Imaging Atmospheric 
Cherenkov Telescope experiments by the Fermi-LAT collabo- 
ration on 27th October 2009 (Fermi-LAT Collaboration 2009, 
priv. comm.). B3 2247+381 is listed in the first Fermi-LAT cat- 
alog of AGNs (Abdo et al. 2010) as 1FGL J2250. 1 +3825 with a 
very hard spectrum (spectral index of -1.6 + 0.1). In the second 
Fermi-LAT catalog this object is listed as 2FGL J2250+3825 
with a spectral index of -1.84 + 0.11 (Ackerman, M., Ajello, M., 
Allafort, A., Antolini, E. et al., ApJ, in prep.). Neronov et al. 
(2010) found a hint of signal at 2.73 cr in the Fermi-LAT data 
above lOOGeV over the period August 2008- April 2010. 

In this paper we present the first detection of VHE y-ray 
emission triggered by the optical high state of the source and 
the first optical light curve of B3 2247+381. We also present si- 
multaneous X-ray data obtained by the Swift and High energy 
(HE) data obtained by Fermi-LAT satellites. 

2. Observations and data analysis 

2.1. MAGIC observations and data analysis 

The VHE y-ray observations were carried out with the MAGIC 
telescopes located on the Canary Island of La Palma (28.8° N, 
17.8° W at 2200 m.a.s.l). The two 17 m telescopes use the atmo- 
spheric Cherenkov imaging technique and allow for measure- 
ments at a threshold as low as 50 GeV 

B3 2247+381 was observed with the MAGIC telescopes dur- 
ing 13 nights between September 30th and October 30th 2010 
collecting a total of 21.2 hours of data, of which 5.3 hours were 
discarded based on the event rate. The effective time of this ob- 
servation, correcting for the dead time of the trigger and readout 
systems is 14.2 hours. Part of the data was taken under moderate 
moonlight conditions. 

All the data were taken in the false-source tracking (wobble) 
mode (Fomin et al. 1994), in which the telescope was alternated 
every 20 minutes between two sky positions at 0.4° offset from 
the source. The zenith angle was between 8° and 35°. 

The data was analysed using the standard MAGIC analy- 
sis framework "MARS" as described in Moralejo et al. (2009) 
with additional adaptations incorporating the stereoscopic ob- 
servations. The images were cleaned using timing information 
as described in Aliu et al. (2009) with absolute cleaning levels 
of 6 phe (so-called "core pixels") and 3 phe ("boundary pixels") 
for the first telescope and 9 phe and 4.5 phe respectively for the 
second telescope. The images were parameterized in each tele- 
scope separately following the prescription of Hillas (1985). 

2 Note, that Veron-Cetty & Veron (2006) list the source as a BL Lac 
candidate and that the classification has not been confirmed. 

3 The stacked dataset of all X-ray selected blazars resulted in a sig- 
nificant y-ray excess, which hints towards the fact that the sources were 
emitting at a very low level. 



In order to reconstruct the shower arrival direction we used 
the random forest regression method (RF DISP method, Aleksic 
et al. 2010) which was extended using stereoscopic information 
such as the height of the shower maximum and the impact dis- 
tance of the shower on the ground. We estimated the RF DISP 
for each telescope separately and obtained two possible solutions 
along the major axis of the shower image in each telescope, re- 
spectively. Finally, we searched for the combination of two so- 
lutions (one from each telescope) that have the shortest squared 
angular distance between them. If this squared angular distance 
is greater than 0.05 degree 2 the event is removed from further 
analysis. This improves the background rejection since hadron 
induced showers have a larger error on the reconstruction of the 
arrival direction. The final arrival direction is the mean of the 
two solutions (weighted by the number of pixels in each shower 
image). 

For the gamma-hadron separation the random forest method 
is used (Albert et al. 2008b). In the stereoscopic analysis im- 
age parameters of both telescopes are used as well as the shower 
impact point and the shower height maximum. A detailed de- 
scription of the stereoscopic MAGIC analysis can be found in 
(Aleksic J., Alvarez, E.A., Antonelli, L.A., Antoranz, P. et al., 
Astr. Particle Phys. in prep.). 

2.2. Optical observations and data analysis 

B3 2247+381 has been observed by the Tuorla groupQ, using 
the 35 cm KVA telescope, located on La Palma, since 2006 
(see Aleksic et al. 2011a for a description of the telescope). 
Observations have been made in the R-band. The brightness of 
the object was measured based on stars in the same CCD-frames 
as B3 2247+38 1 . These stars were calibrated by the Tuorla group 
during the observing seasons. R-band magnitudes were con- 
verted to fluxes, using: F(Jy)= 3080 • 10™ mR ^ 25 . The fluxes were 
corrected for galactic absorption by R=0.398 mag (Schlegel et 
al. 1998). During the years 2006-2009 B3 2247+381 was a quite 
faint and steady source at R ~1.8mJy, but during late summer 
2010 it went to a high optical state, reaching an average flux 
level of 2.4 mJy. The source also stayed at this level throughout 
the observing season (see Figure 3). In late September an alert 
was given to MAGIC and it started observations on September 
30th. An alert is issued when the optical flux has increased by 
50% from the long term running average. 

2.3. Swift observations and data analysis 

The prime objective of the Swift Gamma-Ray Burst observa- 
tory, launched in November 2004 (Gehrels et al. 2004), was to 
detect and follow up on Gamma-Ray bursts, but it has turned 
into a multi-purpose observatory due to its fast slewing and 
response capacity and its multi-wavelength coverage. Swift is 
equipped with three telescopes, the Burst Alert Telescope (BAT; 
Barthelmy 2005), which covers the 15-150 keV range, the X-ray 
telescope (XRT; Burrows et al. 2005) covering the 0.3-10keV 
energy band, and the UV/Optical Telescope (UVOT; Roming et 
al. 2005) covering the 1800-6000 A wavelength range. 

Swift Target of Opportunity observations were requested and 
from October, 5 to 16, 2010. SwiftPiRT observed the source for 
~ 5 ks every night, in Photon counting (PC) mode. We also anal- 
ysed Swift archival data from August 10th, 2009, February 18th, 
2010 and April 18th, 2010 in order to compare the level of the 
X-ray emission to previous observations. 

4 http://users.utu.fi/kani/ 
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The data were processed with standard procedures using the 
FTOOLS task XRTPIPELINE (version 0.12.6) distributed by 
HEASARC within the HEASoft package (v.6.10). Events with 
grades 0-12 were selected for the PC data (see Burrows et al. 
2005) and the response matrices available in the Swift CALDB 
(v.20 100802) were used. For the spectral analysis, we extracted 
the PC source events in the 0.3-10keV range within a circle with 
a radius of 20 pixels (~ 47 arcsec). The background was ex- 
tracted from an off-source circular region of 40-pixels radius. 

The spectra were extracted from the corresponding event 
files and binned using GRPPHA to ensure a minimum of 
30 counts per bin in a manner so that the x 2 statistic could 
reliably be used. Spectral analyses were performed using 
XSPEC version 12.6.0. We adopted both a simple power law 
model and a log-parabolic model as in Massaro et al. (2004) 
with an absorption hydrogen-equivalent column density fixed 
to the Galactic value in the direction of the source, namely 
1 .2 x 10 21 crrT 2 . The two models provide similar results in terms 
of goodness of fit above ~ 0.7 keV. However below this energy 
the differences are in general negligible due to low statistics. 



g 350 F" 



9 distribution 



100 - 



50 - 



329.0± 18.1 



Time = 14.23 h 
N on = 488; N otf . 
N„ = 159.0 

Significance (Li&Ma) = 5.58o 




0.1 



0.2 



0.3 



0.4 
i 2 [degrees 2 ] 



Fig. 1. Distribution of the squared angular distance (0 2 ) for the 
on-source counts in the direction of B3 2247+381 (black points 
with errorbars) and the normalized off-source events (gray his- 
togram). 



Swift/UVOT observed the source in the "filter of the day" 
mode, that is a different filter was used for different observations. 
This does not allow to compare the UV fluxes among different 
days. UVOT source counts were extracted from a circular region 
5 arcsec-sized centered on the source position, while the back- 
ground was extracted from a larger circular nearby sourcefree 
region. These data were processed with the uvotmaghist task 
of the HEASOFT package. The observed magnitudes have been 
corrected for Galactic extinction Eg-v - 0.149 mag (Schlegel et 
al. 1998). 

2.4. Fermi data analysis 

The Fermz-LAT is a pair conversion telescope designed to 
cover the energy band from 20MeV to greater than 300 GeV 
(Atwood et al. 2009) which operates in survey mode, scan- 
ning the entire sky every three hours. The data used in this 
paper encompasses the time interval from August 5th, 2008 
to April 7th, 2011 (MJD 54683 - 55658), and were analyzed 
with the Fermi Science Tools package version v9r23p0, which 
are available from the Fermi Science Support Center (FSSC). 
Only events with the highest probability of being photons, 
those in the diffuse class, located within 12° of B3 2247+304 
were used in this analysis. In addition, a cut on the maxi- 
mum zenith angle (< 100°) was applied to reduce the con- 
tamination from the Earth-limb gamma-rays, which is pro- 
duced by cosmic-rays interacting with the upper atmosphere. 
The background model used to extract the y-ray signal includes 
a Galactic diffuse emission component and an isotropic com- 
ponent (including residual instrument background), modelled 
with the files gllJem_v02_P6.Vll_DIFFUSE.fit and isotropic 
iem_v02_P6_Vll_DIFFUSE.txt, which are publicly available^. 
The normalizations of the components comprising the total 
background model were allowed to vary freely during the spec- 
tral point fitting. The spectral fluxes were derived with the post- 
launch instrument response functions (IRF) P6_V 11 .DIFFUSE, 
and applying an unbinned maximum likelihood technique 
(Mattox et al. 1996) to events in the energy range spanning from 
300 Me V to 300 GeV. All the sources from the 2FGL catalog 
(Abdo, A.A., Ackermann, M., Ajello, M., Allafort, A. et al., 
ApJS, in prep.) located within 7° radius of B3 2247+38 were 



included in the model of the region. The initial parameters in 
the XML file were set to those of the 2FGL catalog, leaving the 
normalization parameters free in the fitting procedure. The sys- 
tematic uncertainty in the flux is estimated as 10% at 100 Me V, 
5% at 560 MeV and 20% at 10 GeV and abovefl 

For the period of the MAGIC observations (30 days between 
September 30th and October 30th 2010), the source is not sig- 
nificantly resolved, and hence only 95% confidence level upper 
limits were obtained. In the light curve presented in Fig. 3, for 
each time bin, if the TS value for the source was TS < 4 or the 
number of predicted photons N-pred <3, the values of the fluxes 
were replaced by the 95% confidence level upper limits. The 2- 
sigma upper limits were computed using the Bayesian method 
(Helene 1983), where the likelihood is integrated from up to 
the flux that encompasses 95% of the posterior probability. 

3. Results 

In the MAGIC data the distribution of the square of the angle 
between the reconstructed direction of the events after cuts and 
the position of B3 2247+381 (RA: 22.83472°, DEC: 38.41028°, 
J2000 as in Ficarra et al. 1985) shows an excess of 159 ± 28 
events above a threshold of 200 GeV. The cuts were previously 
optimized on a trial sample of Crab Nebula data. The threshold 
was calculated from Monte Carlo simulated data by finding the 
peak of the differential rate vs. energy distribution after cuts and 
spectral re-weighting. The expected background level with the 
same cuts is 329 events, calculated from the 8 2 distribution of 
the reconstructed direction of the events with respect to the anti- 
source position, located 180° from the real source position in the 
camera plane (Figure 1). The measured excess corresponds to a 
post-trial significance of 5.6 cr calculated using eq. 17 from Li 
and Ma (1983). 

The source position and extension, determined by a 2D 
Gaussian fit of the sky map produced with the cuts above, are 
consistent with a point-like source placed at the position of B3 
2247+381 within 0.015°, well within the statistical uncertainty 
and the systematic pointing uncertainties of MAGIC. 

The integral flux of the source above 200 GeV was estimated 
to be (5.0 + 0.6 stat + l.l sys ) • 10~ 12 ph cirrV 1 . The differen- 



5 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 6 http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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Fig. 2. The unfolded differential energy spectrum of B3 
2247+381 observed by MAGIC. The black data points refer to 
the measured spectrum, while the grey dashed points have been 
corrected for the attenuation of the EBL. The thick black and 
dashed grey lines are power-law fits to the respective data points 
(fit results are given in the text). The dashed band corresponds 
to the statistical error of the fit to the measured spectrum, while 
the white band surrounding it is the sum of the statistical and 
systematic errors of the fit. 



tial energy spectrum is well described by a simple power-law: 
dN/dE = f (E/300GeVF. The photon index y was found 



to be -3.2 + 0.5 S 



0.5 sys , and a flux normalization ft at 



In 



300GeV of (1.4 ± 0.3 stat ± 0.2 sys ) • 10" 11 ph cnrV^eV" 1 . 
order to correct the effects in the spectrum determination intro- 
duced by the limited energy resolution, different unfolding al- 
gorithms (Forward, Tikhonov, Schmelling 1&2, Bertero; all de- 
scribed in Albert et al. (2007c)) were used, and all agreed within 
errors. Taking into account the attenuation due to pair produc- 
tion with the extragalactic background light (EBL), the spec- 
trum was found to be compatible with a power law with pho- 
ton index y = -2.7 ± 0.5 stat ± 0.5 sys and flux at 300 GeV ft = 
(2.0 + 0.3 smt + 0.3 sys ) ■ 10- 11 ph cnrV'TeV- 1 (Figure 2). Two 
different EBL models were used (Dommguez et al. 2011 and 
Kneiske & Dole 2010), and they were found to be in good agree- 
ment with each other, well within the statistical uncertainties. 

Long term light curves of B3 2247+381 in VHE y-rays 
(MAGIC), HE y-rays (Fermi-LKT), X-rays (Swift) and optical 
(Tuorla Observatory) are shown in Figure 3. Our detection in 
VHE y-rays is compatible with the previous upper limit from 
2006 and thus no variability can be established in this energy 
band. However, in X-rays and the optical band a clear increase 
of the flux in fall 2010 is evident, while the Fermz'-LAT light 
curve is consistent with a constant flux. A fit with a constant to 
the eleven Fermi flux points where the source is significantly de- 
tected, gave a flux value of (3.7 + 0.5) • 10~ 9 ph cm s , with a 
reduced ;f 2 of 0.7 with ten degrees of freedom. The Fermz'-LAT 
is not sensitive enough for detecting short term variations at this 
flux level. The temporal evolution of the VHE y-ray, X-ray and 
optical emission from B3 2247+381 during September-October 
2010 observation shows no strong variability on time scales of 
a night (Figure 4). In particular, the MAGIC light curve above 
200 GeV is consistent with a non-variable source, having a re- 
duced x 2 of 0.6 with ten degrees of freedom. During one night 
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Fig. 3. Long term light curves of B3 2247+381 in VHE y-rays, 
Fermi-LAT HE y-rays (two months time intervals), Swift X-rays 
and optical KVA R-band. 



the X-ray flux is significantly higher (almost factor 2) than the 
other X-ray points, but unfortunately we do not have simultane- 
ous optical or MAGIC data for that night. 



4. Discussion 

In this paper we report the discovery of VHE y-rays from 
B3 2247+381 by MAGIC. The MAGIC observations were trig- 
gered by an optical high state of the source, like several other 
discoveries. However, the observed VHE y-ray flux is consistent 
with the upper limit from previous observations and we there- 
fore cannot conclude if the source was in a higher VHE y-ray 
state during the observations. 

In Figure 5 we show the spectral energy distribution of the 
source during the MAGIC observations, together with simulta- 
neous Swift and optical data, and other non-simultaneous data. 
The Swift observations show that the source was in a high state 
also in X-rays. In the Fermi energy range the source is very 
weak, which limits the capability of detecting statistically sig- 
nificant flux-variability on time scales of a few months. The syn- 
chrotron component of the SED is showing a significantly larger 
emission in the high state, while the inverse Compton compo- 
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Fig. 4. Same as Fig. 3, but zoomed into the time interval of the 
MAGIC observations in September-October 2010. 
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Fig. 5. Spectral energy distribution of B3 2247+381 (red: EBL 
corrected MAGIC spectral points). The green crosses are 1FGL 
Fermi data points (Abdo et al. 2010), while the pink points rep- 
resent the Fermi analysis from this work (2.5 years of data). 
Blue arrows show the 95% confidence upper limits computed 
from Fermi-LAT data for the time interval of the MAGIC ob- 
servation. Low (high) state Swift data were taken on April 18th 
2010 (October 5-16, 2010). The host galaxy contribution has 
been subtracted from the KVA R-band data (red and light blue 
squares), following Nilsson et al. (2007). The data have been 
corrected for galactic absorption. Green and light blue points 
represent non-simultaneous low state data. The solid line is our 
SSC-model fit to the high state observations; the dotted line is a 
fit to the low state observations. 



nent is consistent with only minor changes. We must however 
note that the weak detection in the HE and VHE y-ray band sig- 
nificantly limits the determination of the inverse Compton com- 
ponent. 

We reproduce the SED with a one-zone synchrotron-self- 
Compton (SSC) model (see Tavecchio et al. 2001 for a descrip- 
tion). In brief, the emission region is assumed to be spherical, 
with a radius R, filled with a tangled magnetic field of intensity 
B. The relativistic electrons follow a smoothed broken power- 
law energy distribution specified by the limits y m i n , y max and the 
break at y, as well as the slopes n\ and tii before and after the 
break, respectively. Relativistic effects are taken into account by 
the Doppler factor 5. The used input model parameters are shown 
in Table 1 . 

To reproduce the change of the Compton and synchrotron lu- 
minosity ratio between the low and the high state we mainly act 
on the electron normalization, source radius and Doppler factor 
(with slight changes also to the other parameters). The steeper 
X-ray slope in the low state implies a larger value of «2- 

The comparison with parameters derived for BL Lac objects 
(e.g. Tavecchio et al. 2010) reveals that the parameters used 
for B3 2247+381 are close to the typical values. As for other 
sources, the somewhat larger (lower) value of the Doppler factor 



6 (the magnetic field intensity B) with respect to "standard" val- 
ues is mainly due to the relatively large separation between the 
synchrotron and IC peaks. 

Additionally, we model the SED, during the high state, with 
the one-zone SSC code from Weidinger et al. (2010) which is 
shown in Figure 6. In contrast to the model from Tavecchio et 
al. (2001) all parameters are basic physical parameters and the 
electron and photon spectrum and their breaks are derived self- 
consistently with a continuous injection of monochromatic elec- 
trons at the energy yo = 10 4 and injection rate K = 8.4 • 10 4 
crrT 3 s _1 . The spectrum is the resulting steady-state solution. 
The environment is defined by the magnetic field B = 0.07 G, 
the acceleration efficiency t acc /t esc = 1.09 (i.e. the particle spec- 
tral index is s = 2.09 and the resulting y max = 4.8 • 10 5 ), and 
the blob radius R — 1.3 • 10 16 cm. The break in the e~ spectrum 
of 1 at jb = 2.9 • 10 4 arises self-consistently from IC and syn- 
chrotron cooling. The common parameters of both models agree 
very well. 

The optical monitoring of candidate VHE y-ray blazars has 
proven to be a successful tool for their discovery. However, for 
B3 2247+381 (like for Mrk 180 and 1ES 1011+496) we can- 
not establish a connection between the optical high state and the 
VHE y-ray emission since the upper limit from previous obser- 
vations (during a low optical state) is higher than the detected 
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Table 1. Input parameters for the high and low states of the SSC-model shown as solid and dashed lines in Fig. 5. For more 
explanations see text. 



Flux State 




Tb 




n\ 


n 2 


B 


K 


S 


R 














G 






cm 


High 


3 -10 j 


7.1 • 10 4 


6 -HP 


2.0 


4.35 


0.06 


2.5 -10 3 


35 


8- lO 13 


Low 


3 -10 3 


6.8 • 10 4 


5 • 10 5 


2.0 


5.35 


0.08 


1.15 ■ 10 4 


30 


4- 10 15 




KVA 

Swift XRT 

MAGIC (de-absorbed) 
Fernii LAT 

model SED (high-state) 



10 10 10" 

v(Hz) 



Fig. 6. The solid line shows the SED model using Weidinger et 
al. (2010). The data points are described in the inlay of the figure. 
The fit parameters can be found in the text. 



VHE y-ray flux during the discovery. Therefore further observa- 
tions are still needed to study the connection between these two 
wavebands (see also Reinthal et al. 201 1). 
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